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Abstract: We consider a thermalization process in a 2-dimensional CFT that has a 
holographic description in terms of the gravitational collapse of a thin shell of null dust. 
This model represents a sudden perturbation of the CFT vacuum that communicates a 
uniform energy density to the system. We study the evolution of two-point functions at 
spacelike separated points (t\, I) and fa, 0), and reproduce the generic pattern first derived 
from the analysis of quantum quenches to critical systems. A crucial characteristic of these 
setups is that the excitations generated by the initial perturbation presents non-trivial 
quantum correlations. As a consequence, for any t j < oo equilibration is only effective on 
finite regions whose size grows as a lightfront. The behavior on larger regions is greatly 
determined by the initial state, which for the quenches we consider and the holographic 
model has relevant differences. However in both cases for late times the dependence on 
the scale I of the two-point functions enters through the effective distance l — t\ — t<2.- We 
interpret the onset of this behavior as an equilibration time for occupation numbers in 
these 2-dimensional models. 
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1. Introduction 

The AdS/CFT correspondence represents a conceptual breakthrough in the study of strongly 
coupled field theories. It has proven to be extremely valuable to obtain universal properties 
of field theories in the strongly coupled regime, although the concrete models accessible 
to construction are often not as realistic as desirable. The correspondence has been very 
successfully applied in recent years to the study of the transport properties of strongly 
coupled plasmas in the context of linear [1] and non- linear [2,3] fluid dynamics. A natural 
continuation of this line of work, which is attracting a growing interest, is its application 
to the study of thermalization processes starting from a far from equilibrium state [4-17]. 

The holographic dictionary relates the plasma phase of strongly coupled field theories 
at thermal equilibrium to black hole geometries [18]. Hence the holographic counterpart 
of a thermalization process must correspond to a process of gravitational collapse ending 
in the formation of a black hole [19-22]. The simplest observables which can be used as 
probes of the thermalization process are one-point functions, since they can be derived 
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from an expansion of the dual background close to the boundary. In [9, 10, 12, 15] it was 
initiated the study of the holographic evolution using observables which need information 
from the dual geometry far from the boundary. 

The dual geometry used in [9, 10, 12, 15] was a Vaidya metric. These metrics describe 
the collapse of a shell of null dust to form a black hole, and are known analytically also 
for spacetimes with negative cosmological constant. This setup represents a CFT starting 
in its vacuum state and undergoing a translationally invariant perturbation, modeled by 
the collapsing shell, which brings it out of equilibrium. After the perturbation ceases, the 
system evolves according to the initial CFT hamiltonian. A similar dynamical setup is 
provided by a quantum quench, which denotes an action on a system in which a parameter 
of the hamiltonian is suddenly changed triggering the subsequent evolution. Holographic 
models for quantum quenches that affect some localized degrees of freedom have been 
proposed in [8, 11]. 

Quantum quenches are unitary processes. As such, the excitations produced by the 
quench present non-trivial quantum correlations. This has very important consequences for 
the evolution towards equilibration of the system. In order that an observable involving 
an scale I has reached its equilibrium value at a time t, it must happen that quantum 
entangled excitations have separated at least a distance I at that time [23] . Causality then 
implies that relaxation does not happen in the system globally, but it takes longer the 
bigger the region considered. Thermalization processes described holographically through 
gravitational collapse are also unitary [9,24]. Using a 3-dimensional Vaidya metric, the 
dependence of the size of thermalized regions with time was reproduced in [9] in perfect 
agreement with previous results [23] . The observable used to follow the dynamical process 
in both works was the entanglement entropy. 

There are some relevant differences between the Vaidya model of thermalization and 
quantum quenches. Quenches from a gapped to a critical system were analyzed at a 
general level in [23,25,26]. This type of quenches differ from the holographic model in the 
entanglement pattern of the initial state that triggers the evolution. Namely in the former 
entanglement is initially localized on small scales, whereas in the Vaidya model there are 
long range correlations in the early stages of the evolution [9] . A quantum quench between 
two different critical models was studied in [27], exhibiting power law correlations at late 
times after the quench. The holographic model has in common with it the existence of long 
range correlations in the initial state. However at late times it presents by construction a 
thermal behavior, in closer analogy with quenches from gapped models. 

The following remark is in order. Although in a unitary process information is not lost 
at the microscopical level, realistic perturbations in extended systems generically populate 
many energy levels and lead to relaxation towards an state that at the macroscopical level 
can be described as thermal [28]. This is the case of the holographic models based on 
gravitational collapse. The final state towards which a system relaxes can however differ 
from thermal equilibrium if additional integrals of motion are present. This happens for 
example in the quench between two critical systems of [27], and in the quench from a 
massive to a massless free boson treated in [26]. 

In [10] the holographic evolution of entanglement entropy was considered for boundary 
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dimension d = 3, with analogous results to [9]. In addition to the entanglement entropy, 
equal-time two point functions and Wilson loops were studied in [12, 15] as probes of a 
thermalization process for d = 2,3, 4. The dynamical background they considered was the 
collapse of an infinitely thin shell of null dust, a limit of the Vaidya setup. For d = 2, 
this simplified background allowed to obtain exact expressions instead of the numerical 
methods necessary in the smooth Vaidya case. 

The aim of this paper is to calculate general two-point functions in the thin shell geom- 
etry. We will compare our results with those of [25,26], where the evolution of one and two- 
point functions for quenches from gapped to gapless systems was analyzed, finding again 
perfect agreement. The general methods used in [23,25,26] to study the dynamics of these 
quenches were based on the special properties of the conformal group in two-dimensions. 
They rely in performing the calculations on euclidean signature, and then prolonging to 
real time. The holographic methods allow for a direct calculation on lorentzian signature. 
Although we are using a very simple geometry to simulate a thermalization process, it is 
likely that our conclusions apply equally to more general models of thermalization based 
on gravitational collapse. 

Besides the influence of the initial pattern of quantum entanglement, a very important 
aspect of the evolution towards equilibration is the redistribution of energy among the 
excitations generated by the initial perturbation. It was suggested in [9] the possibility of 
deriving from the holographic evolution of the entanglement entropy a time at which the 
redistribution of energy among modes according to thermal equilibrium should be nearly 
completed. An important result of the present work is to show that precisely the same 
threshold time is obtained from the analysis of general two-point functions. 

The organization of the paper is as follows. Section 2 summarizes the results for the 
evolution of two-point functions after a quantum quench derived in [25,26]. In Section 3 we 
recall the holographic evaluation of two-point functions for operators with high conformal 
dimension, that can be obtained in terms of the proper length of geodesies that anchor 
in the AdS boundary [19,32]. We review the construction of geodesies in AdS3 and in a 
BTZ black hole. The derivation of geodesies in the infalling shell geometry is addressed in 
Section 4. Sections 5 and 6 consider particular limits in which explicit expressions for the 
holographic two-point functions can be obtained: the case of large space separations for 
the two-point functions and that of late times after the perturbation respectively. Section 7 
contains a discussion of our results. Several technical facts are collected in two appendices. 



2. Two-point functions after a quantum quench 

The two-point function of primary operators in the vacuum of an euclidean 2-dimensional 
CFT on the plane is completely fixed by symmetries. Since we will be interested in dy- 
namical situations we need its prolongation to real time, with the result 

(O(h,l)O(t 2 ,0)) = \ (2.1) 
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where A is the conformal dimension of the operator O and At = ii — t 2 - In the same way, 
conformal symmetry determines the two-point function at thermal equilibrium 



(0(Mi)0(O,i2)) 



2ttV 



P) 2(cosh(2vr///?) -cosh(2vrAi/^)) 



(2.2) 



where (3 is the inverse of the temperature. Both in vacuum and at thermal equilibrium the 
two-point functions have singularities on the light-cone, meaning that in a 2-dimensional 
CFT all excitations propagate at the speed of light even in the presence of a thermal 
bath. In order to make precise the prolongation to real time of the euclidean propagator, 
a prescription has to be given that determines how to treat the lightcone singularities. 
Depending on the prescription chosen, the Wightman function or the time-ordered two- 
point function can be obtained. The Wightman functions above require At—>A.t — i0 + . 

The evolution of one- and two-point functions after a quantum quench was analyzed 
in [25,26], following previous work on the entanglement entropy [23]. As in that case, 
explicit results can be derived for 2-dimensional field theories when the system before the 
quench is in the ground state of a massive hamiltonian, while the dynamics after the quench 
is conformal. The initial ground state, which is not an eigenstate of the new hamiltonian, 
provides the starting point for the evolution. As customary, the quench is taken to happen 
at t = 0. 

In [23,25,26] it was argued that the evolution after a quantum quench can be obtained 
from the continuation to real time of two-point functions in a strip geometry of broadness 
2to, with ro a parameter of the order of the inverse mass gap in the initial state. The strip 
geometry can be mapped to the upper half-plane through the conformal transformation 
oj(z) = ^ log z, where cj and z are the coordinates in the strip and upper half-plane (UHP) 
respectively. Then 

(O(n,/)O(r 2 ,0)) = \uj'(z 1 )\- a \uj'(z 2 )\- a (0(z 1 )0(z 2 )) vrp , (2.3) 
with uj(zi) = I +*Ti and uj(z 2 ) = ir 2 . The two-point function on UHP is well known [29] 



(- 
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<0(3l)0(* 2 )>UHP = ( 

\Z\2Z\ 2 Z X 

where Zij = \z{ — Zj\ and the function F depends on the particular CFT and boundary 
conditions under consideration. Rotating to real time, Tj = To+itj, the model independent 
piece of the two-point function fl2.3|) is 



(2.5) 



7r \ 2 cosh(7rZ/2ro) + cosh(7r(ti + t2)/2ro) 

2tq J 4 cosh (-7rti /2ro ) cosh (7rt2/2ro) (cosh (-7r//2ro) — cosh (71- (ti — £2)/2to)) 

Taking moreover the limit ti, 1^>tq, the argument of the function F gives 

Zll Z 2 2 e ^l+*2)/2r 



12^21 



e irl/2To _|_ e 7r(ti+t2)/2ro 



(2.6) 
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When I < ti+t 2 the argument of F is approximately one, and for any CFT it is verified 
F(l) = 1 [25,26]. In this case, the contribution ( |2.5| ) reproduces the result for the two- 
point function at thermal equilibrium ( |2,2j ) with inverse temperature j3 = 4tq. This leads 
to interpret 1/4to as an effective temperature for the long wavelength modes. This same 
conclusion was obtained before from the evolution of the entanglement entropy [23]. 

When instead I > t±+t 2 , the argument of F is very small, of order e ~ 7r (' - *i - *2)/ 2 ' r o. 
Remarkably, also in this limit F has a generic form [25,26,29] 

F(z) oc z x . (2.7) 

The exponent x depends both on the operator O and the initial state for the evolution, 
and it is only different from zero if O has a vanishing one-point function. In that case, and 
assuming again large times and interval we have 

(0(1, h)O(0, t 2 )) OC e -A7r(ti+t 2 )/2ro e - :r 7r(/-ti-t2)/2To _ ^.8) 

The transition region between both regimes happens at l^t\+t 2 over a region of order To- 
It is governed by the function F and thus model dependent. 

The physical picture that explains the previous evolution pattern is as follows. The 
non-zero energy density created at t = by the quench translates into the production of 
a translationally invariant sea of excitations. As the quench is a unitary process, there 
will be quantum correlation between the resulting excitations. Entangled excitations will 
tend to spread over ever larger regions. In order that observables inside a region can 
behave as thermalized, the distance between left and right moving entangled excitations 
must be bigger than the size of the region. In this way each component of an entangled 
pair can act as part of a thermal bath for the other. A simple picture where entangled 
excitations are assumed to move classically and without scattering, was further proposed 
in [23]. For the critical quench described above, all excitations are taken to propagate at 
the speed of light. Since the ground state before the quench has a mass gap and hence a 
finite correlation length, only excitations produced approximately at the same point will 
be entangled. Both facts give rise to a sharp horizon effect. In order that a region of size 
I behaves as thermalized, a time 1/2 after the quench is required. Moreover a variation in 
the size of a region of the order of the initial correlation length around I = 2t will bring 
from regions that behave as thermalized to those that contain the main information about 
the entanglement pattern of the initial state. In this second case, the expectation value of 
observables should be reminiscent of that in the initial state. This is indeed what happens 
for the entanglement entropy, which is defined as the von Neumann entropy of the reduced 
density matrix associated to a subregion of the system. Its evolution for the case of a single 
interval was studied in [23]. 

The same picture explains the behavior of non-equal time two-point functions [25,26]. 
The distance that separates a pair of left and right moving entangled excitations produced 
at t = and propagating up to generic different times is t\ + t 2 . When I < t\ + t 2 there 
are no entanglement excitations over the distance involved by the two-point function and 
a thermal result is reproduced. On the contrary, when I > ti + t 2 there are entangled 
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excitations on distances of order I. Their pattern of entanglement is determined by the 
initial state, while the effective distance their quantum correlations feel is l — t\—t2- The 
fact that their contribution to the two-point function ( |2.8| ) is exponentially suppressed 
can be traced back to the fact that the initial state is massive. Notice that x/tq ~ xm, 
where m is the mass gap in the initial state and a; is a parameter containing information 
about the initial state. The two-point function ( |2.8| ) has also an exponential suppression 
with t± + t2- This combination gives a measure of the decoherence due to the growing 
separation of entangled excitations. This fact is independent of the initial state. Hence the 
only information about the initial state involved in it is to, which sets the energy density 
generated by the quench. 

3. Holographic two-point functions 

According to the AdS/CFT dictionary, two-point functions of gauge invariant operators are 
derived from the on-shell supergravity action for the bulk fields that couple to the chosen 
operators [30,31]. This requires to know the bulk to boundary propagator associated to 
that field, which can only be obtained analytically in very simple examples. In particular, 
for a time-dependent background the derivation of the bulk to boundary propagator will 
involve solving numerically a partial differential equation. 

An alternative way to derive boundary Green's functions was proposed in [19]. They 
should be given in terms of bulk correlators for the associated fields to which the operators 
couple. Bulk propagators can be represented as an integral over paths that join the insertion 
points (ii,Z;r&) and (i2, 0; n,), having then 

(O(t u l)O(t 2 ,0)) = lim r 2 b A [ VVe~ AL{v \ (3.1) 

where L(V) is the proper length of the path with the convention to be real for space-like 
trajectories. For operators with very high conformal dimension A it is justified to perform 
a saddle point approximation, which reduces the path integral to a sum over geodesies [32]. 
It should be stressed that this approximation requires special care in lorentzian spaces 
because of subtleties associated with the choice of steepest descendent contour [33-35] . 
We will assume that these complications do not arise in the 3-dimensional collapsing shell 
background we will be interested in, as it was also done in [12, 15]. The consistency of the 
results we will serve as a check of this assumption. 

The geodesic length has a divergent logarithmic contribution from the region r b — >oo. 
A regularized geodesic length can be defined by subtracting this piece 

L = L reg + 2logr b . (3.2) 

The logarithmic contribution to L precisely cancels the power of rf, in ( |3.1| ) and the limit in 
this expression can be easily taken. Relation ( |3.1j ) reduces then to evaluating the regularized 

1 It was shown in [34] that for Schwarzchild-AdSd, with d > 3, there exist 2-point functions receiving 
contributions from several geodesies. In some cases the largest contribution comes not from a real but from 
a complexified geodesic. This analysis requires the geometry of the space under consideration to admit an 
analytic continuation, which is not the case of the infinitely thin shell background. 
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length of geodesies with endpoints in the boundary. In the 3-dimensional case we are 
considering, there is a unique geodesic that joins two given boundary points. Hence the 
two-point function of operators with high conformal dimension will be simply given by 

(0(ti,l)0(t 2 , 0)> = e -AW*i.i;fe,o) _ (3_ 3) 

Two-point functions are not the only observable whose holographic determination in- 
volves evaluating the length of bulk geodesies. This is also the case of the entanglement 
entropy [36,37]. The entanglement entropy of a single interval in a 2-dimensional CFT is 
given holographically by the length of the equal-time geodesic that anchors on the bound- 
ary at the interval endpoints. The study of equal-time geodesies in a Vaidya geometry 
representing the collapse of a null dust shell with finite thickness was carried out numeri- 
cally in [9]. The limit of an infinitely thin shell in AdS3 allows to use analytic instead of 
numerical methods. Equal-time geodesies were reconsidered in [12, 15] in this background. 
We want to extend this study to general geodesies with endpoints in the boundary. Us- 
ing the approximation (|3.3| ), we will reproduce the pattern for the evolution of two-point 
functions presented in the previous section. 

The geometry describing the infinitely thin shell is that of a BTZ black hole outside 
the shell and AdS3 inside. Both spaces join along a radial null surface that we take to 
start at the boundary at t = in analogy with the convention for the quantum quench. 
This geometry represents an instantaneous action on the dual CFT vacuum that brings 
the system out of equilibrium by communicating an homogeneous energy density. The 
system will evolve subsequently towards equilibration at a temperature set by that of the 
dual black hole formed in the gravitational collapse process. In this set up, the only field 
theory observable with a non-zero one-point function is the stress tensor. The operators O 
for which we will compute the two-point functions using the geodesic approximation have 
thus vanishing one-point functions, as it was assumed in the quantum quench result fl2.g|). 

We will review now the description of geodesies in the AdS3 and BTZ geometries, since 
we will need them to evaluate two-point functions in the thin shell background. The metric 
of the non-compact BTZ black hole is given by 

ds 2 = - (r 2 - m) dt 2 + — + r 2 dx 2 , (3.4) 

r A — m 

where the parameter m represents the black hole mass. The associated Hawking tem- 
perature is T = ^/m/2ir and the radius of the event horizon = \pm. The case m = 
of the previous metric describes the AdS3 solution in the Poincare patch. In the affine 
parameterization, the geodesic equations are 

t = (3.5) 

f = ±- v / A 2 r 2 + (r 2 - m) {Br 2 - rj) (3.6) 

x = ^ (3.7) 

where r*, A are integrals of motion associated with the Killing vectors d x and dt, and there- 
fore having an interpretation of conserved momentum and energy. The previous equations 
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imply x^x v gav = B. If we take A to represent the geodesic proper length we have B = 1, 0, — 1 
for spacelike, lightlike and timelike geodesies respectively. Equation ([O]) implies that only 
when B = 1 the associated geodesies can reach the boundary. Since this is what we need 
for the holographic calculation of the CFT two-point functions, from now on we will focus 
on spacelike geodesies in the bulk. 

It will be convenient to describe the geodesies as functions t(r),x(r) of the radial 
coordinate, instead of using the affine parameterization. As it is clear from ( |3.6| ) each 
function will have two branches, which we will denote as ±. We choose them to be defined 
by the conditions 

t+{°°) > £-(oo) , x + (oo) > x_(oo) . (3-8) 

The integration constant r* and A can have either sign. The reason for the criterium ( |3.8[ ) is 
that the resulting functions t±(r) and x±(r) depend only on the modulus of the integration 
constants. Hence without loss of generality we can consider r* and A to be always positive, 
and shift the information about their sign to the choice of branch. Since the signs of both 
constants are independent, so will it be the choice of branch for the functions t and x. 

3.1 AdS 3 

The geodesic equations (p.5[)-( |3T7[ ) are straightforward to integrate for AdS3. The result 
for spacelike geodesies is 

t±(r) = t ± t AdS (r) , x±(r) = x ± x AdS {r) , (3.9) 

where 

tAdS = XAdS = \/A 2 + r 2 -r 2 
A r* r\r 2 - A 2 \ 

Given the translational invariance in t and x of the geometry, the integration constants 
to and xq do not have any influence on the geodesic length. The separation between the 
geodesic endpoints {t\,x\) and (£21^2), is determined by the constants of motion r* and A 
as 

where At = t\— ti and l = \x\— x-z\- 

When r* > A, the geodesies extend from the boundary to the value of the radial 
coordinate 

r m = V / ^42 5 (3-12) 



where both branches join. As can be seen from ( 3.11 ), the geodesic endpoints are spacelike 
separated in this case. The lightlike case can be reached from the spacelike one in the limit 
r* —7- A Keeping the separation between the geodesic endpoints ( |3.11 ) finite in this limit 
requires 

n = p, A = p-l/l, (3.13) 

with p— > 00. The minimal value of the radial coordinate r m tends then to infinity. This 
shows that geodesies joining two lightlike separated points in the case of a non-compact 
boundary, are lightlike geodesies that do not leave the boundary. 
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When instead r* < A the branches ± describe disconnected trajectories that reach 
down to r = 0, where t, x = ±00. In this case there is no reason to take the integration 
constants to an d xo in (|3.9[ ) equal for both branches. Let us however do it, considering 
both trajectories as part of a single geodesic. The separation between its endpoints is again 
given by ( 3.1 1| ) , with the result that they are timelike separated. 



The proper length of the geodesic is 



I -I ? = 2/ , „, 2 y (3.14) 



jro y A 2 + r 2 — r- 



with the lower integration limit equal to r m in ( |3,12 ) for r*> A and zero otherwise. The 



upper integration limit rj, has been introduced in order to avoid the logarithmic divergence 
from the integration region close to the boundary, as explained above. The result for the 
regularized geodesic length (|3,2|), both in the spacelike and timelike cases, is 



r. 



2 



A 2 \ 



L reg = -log^ 1 ^ = log|/ 2 -At 2 |. (3.15) 



Using the geodesic approximation fl3.3j), the vacuum two-point function of primary opera- 
tors with spacelike separation is reproduced [32]. 

It is tempting to apply the geodesic approximation to two-point functions with timelike 
separations. The length of the disconnected geodesies proposed above gives almost the right 
result. Substituting in fl3.3|) , it reproduces the correct value for the two-point function up 
to a constant phase. This mismatch comes from the modulus in the argument of the 
logarithm, which for timelike separations is relevant. The presence of the modulus is 
actually unavoidable because with the convention we have chosen the proper length of 
spacelike geodesies is always a real number. It would be very interesting to understand 
whether the geodesic approximation to the holographic derivation of two-point functions 
can exactly reproduce the timelike case, explaining the mismatch in a phase that we have 
found. It would be also important to determine whether the geodesic approximation can 
represent the different prescriptions for treating the lightcone singularities, leading to the 
Wightman or time-ordered two-point functions. 

3.2 BTZ black hole 

The geodesic equations fl3.5|) -( j3"77l) in a BTZ background can be also explicitly integrated, 
obtaining 



m 1 (eX t (r ) =p sign(r 2 - rn) \J X%{r) - C^j 



t±(r) =t - — log 
x±(r) = x - — log mT 1 (tX x (r) =F V x U r ) ~ c 



(3.16) 



where we have defined 



v , \ 2 a2 2A 2 m 2 a2 2r 2 m 
X t {r)=rl-A 2 -m--„ , X x (r) = r 2 - A + m %- 



(3.17) 



C = (r 2 - A 2 - rn) 2 - AA 2 m = (r 2 - A 2 + m) 2 - Ar 2 m , 
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together with 

e = sign(r 2 — A 2 — m) , e = sign(r 2 — ^4 2 + m) . (3.18) 

The value of the radial coordinate ranges from infinity to the zero of the square root 
in the functions ( |3.16| ) with larger r 2 

r 2 m = ^{r 2 -A 2 + m + Vc), (3.19) 

whenever this quantity is real and positive. Otherwise we have two separated trajectories, 
both reaching the black hole singularity r = 0. A necessary condition for this not to happen 
is that C is positive, which we will assume in this subsection. This is the only case relevant 
for the geodesic approximation to the holographic derivation of two-point functions at 
thermal equilibrium. 

In order to describe the geodesies we will need to consider the extended Penrose dia- 
gram of the BTZ black hole, which has a second asymptotic region (wedge III) and a past 
singularity (wedge IV) . The time t acquires a non- vanishing imaginary part in three of the 
four wedges that compose the extended diagram, as shown in Fig.| [38,39]. 




Figure 1: On the left we show the imaginary part of t(r) in the different wedges of the Penrose 
diagram of the BTZ black hole. On the right we draw qualitatively different types of BTZ geodesies 
according to our classification. 



The different geodesic types can be classified by the signs e and e ( 3.18] ), see also [15]. 



When e = e = 1 the minimal value of the radial coordinate satisfies r m > r^, implying that 
these geodesies stay always outside the event horizon and are contained in wedge I or III. 
When e = — 1 and e = 1 the associated geodesies reach behind the event horizon but they 
do not fall into the singularity, since for them < r m < r^. Close to the horizon the two 
branches of the geodesic time coordinate behave as 

i ± (r)~ T — log(r-r h ). (3.20) 



2 As it should be, the same value is derived from t(r) and x(r). 
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If we choose the branch + to start at the boundary of wedge I, the branch — must start 
at the boundary of wedge III in order that they can join at r m in wedge II 3 . Hence these 
geodesies connect the two separated AdS boundary regions of the eternal black hole (see 
Fig.|l|). Finally when e = e = — 1, we have two disjoint trajectories each of them reaching 
the singularity. As we did for AdS, let us however consider them formally as part of a 
single geodesic characterized by a common value of the integration constants to an d xq. If 
branch + starts on the boundary of wedge I and falls to the future singularity, then branch 
— represents a trajectory that exits the past singularity on wedge IV and returns to wedge 
I. 

The conserved quantities A and r* determine the separation between the endpoints for 
the three types of geodesies described above 

2vrAt \rl-A 2 -m\ 2ttZ \r 2 - A 2 + m\ 
cosh — — = = , cosh — — = = , (3.21) 

P Vc P Vc 

where for geodesies ending on opposite boundaries we have defined At=Re(ii— i2). When 
the endpoints lie on the same boundary their separation is spacelike (timelike) when e = 
e = 1(— 1). For geodesies connecting both AdS boundaries, there is no restriction on the 
values of At and I. 

The proper length of the geodesic in the BTZ background is 

r°° dr f r b dr 2 
L = 2 —= - . (3.22) 

Jr m r J r i m ^ A 2 r 2 + (r 2 - m)(r 2 - r 2 ) 

The integration can be explicitly performed, with the result for the regularized length 



L reg = - log — p . (3.23) 



For geodesies with e = e = — 1, ( 3.19 ) is negative. However the integral in the second 
equality of ( gjgp 

can be generalized to this case by analytically prolonging its integrand 
to negative values of the radial coordinate square r 2 >r^- Remarkably, this leads to the 
same result ( |3.23| ) for the geodesic length. 

Using ffOlD to express C in terms of I and At, and substituting in ( |3.3|) we obtain 



(O(Mi)O(O, h)) 



/2tt 



- 



\P J 2|cosh(2vr///3) - eecosh(2vrAt//3)| 



(3.24) 



When e = e = l, ( 3.24 ) reproduces the two-point function at thermal equilibrium for space- 
like separations in a 2-dimensional CFT [33]. As it was the case in AdS3, e = e = — 1 gives 
the right answer for two-point functions with a timelike separation up to a phase. 

It has been proposed that geodesies connecting the two asymptotic AdS boundaries of 
the eternal black hole are associated with field theory correlators in the Schwinger-Keldysh, 
or real-time formalism [38,40]. In this formalism the tensor product of two copies of the 



3 As already noticed, the integration constants do not need to coincide for both branches. In this case 
we should choose to- =to+—i/3/2, with to+ real. 
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original field theory is considered. The following pure but entangled state is associated 
with the system at thermal equilibrium 

i 

with \Ei) energy eigenstates of the theory. Ordinary thermal correlators are obtained when 
all operators are inserted on the same copy of the system. In the state ( |3.25| ), correlators 
with operators acting on different copies of the system can be related with ordinary ones 
by 

(¥|0i(i,t a )0 2 ((M a )|tf>) = (*|Oi(i,ti)Oi(0,-t 2 -i/9/2)|*). (3.26) 

The holographic 2-point function (|3.24j ) agrees with this relation. Using the bulk to bound- 
ary propagator, two-point function between operators inserted on both copies have been 
obtained in the context of the AdS/CFT correspondence in [38,39]. 



4. Geodesies in the infalling shell geometry 

Our aim is to construct general geodesies with endpoints on the boundary in the dynamical 
geometry describing the collapse of an infinitely thin shell of null dust. In order to describe 
both the metric and the geodesies it is convenient to substitute the time coordinate t, which 
is not constant across the shell, by an infalling radial null coordinate v which it is. Since 
we are taking the shell to start its collapse at t = from the boundary, its trajectory will 
be given by v = 0. Hence the associated metric is 

ds 2 = -{r 2 - m{v))dv 2 + 2drdv + r 2 dx 2 , (4.1) 

with m(v) = for v<0, describing empty AdS3, and m(v) = m for v>0, describing a BTZ 
black hole. The more general Vaidya metrics correspond to consider an arbitrary mass 
function m(v). The relation between v and t in the AdS and BTZ backgrounds is 

v = t-l/r, AdS, (4.2) 

v = t+^-log^^, BTZ. (4.3) 
2r h r + r h 

We will assume that the geodesies endpoints (ti,xi) and (£2)^2) satisfy 

h > t 2 , xi - x 2 = I > 0. (4.4) 

This choice does not represent a loss of generality since parity, which is a symmetry of our 
set up, changes ( f4,4| ) to xi<x 2 - From the criterium (|3.8j), the previous choice implies 

v + (oo) = t\ , V-{oo) = t% , x + (oo) = x\ , x_(oo) = X2 ■ (4-5) 

There are two types of geodesies to consider. We have those that stay entirely outside 
the shell and only feel the BTZ part of the geometry. These geodesies reproduce the result 
for two-point functions at thermal equilibrium. Then there are geodesies that intersect 
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the trajectory of the shell. Actually no matter how large t\ and t<i are taken after the 
perturbation in the case with a non-compact boundary, there will always be an I such 
that the associated geodesic reaches the shell [9]. Indeed the larger I is, the deeper the 
associated geodesies will get in the black hole geometry and in particular the closer to the 
horizon r^. The second term in the rhs of fl4.3| ) makes then a dominant contribution that 
brings the null coordinate v towards small values, untill v = is attained. Geodesies with 
larger values of I will intersect the trajectory of the shell and lead to a non-equilibrium 
result for the two-point functions. 

4.1 Thermal geodesies 

The condition for geodesies to lie in the BTZ part of the geometry was studied in [9] for 
the equal time case t\ = t2 = t. The following very simple relation turns out to hold for 
equal time BTZ geodesies that anchor on the same AdS boundary 

l = 2(t- v m ) , (4.6) 

where v* denotes the minimal value of the coordinate v attained by the geodesic. Since the 
shell follows the trajectory v = 0, geodesies with l<2t will lead to a thermal result for the 
associated observables: entanglement entropy and equal-time two-point functions. This 
reproduces the field theory analysis of quantum quenches reviewed in Section 2. Namely, 
in order that a region of size I can appear as thermalized, a time at least equal to 1/2 has 
to be waited [23] . 

As we also saw in Section 2, the threshold for thermal behavior of two-point functions 
involving different times after a quench was l = t\+t2 with ij > 0. Thus at threshold the 
separation between the insertion points of the two-point function is spacelike. In order to 
reproduce this bound holographically it is then enough to focus on geodesies with spacelike 
separated endpoints, i.e. those contained in wedge I. Equal-time geodesies are characterized 
by having a vanishing energy parameter A, such that the functions v±(r) for them coincide. 
The minimal value of the radial coordinate that they can attain coincides with r* and by 
symmetry at this point t>* is reached. For geodesies with endpoints at different times 
the two branches v±(r) are different, joining smoothly at the minimal value of the radial 
coordinate r m ( 3.19j ), which now differs from r* (see Fig ||). Given the criterium 



chosen to distinguish between branches, it is clear that the minimum value of the null 
coordinate must occur in the branch — . Thus we should search for a zero of 

*L = _JL_fi ^ ^ (47) 

dr r 2 -my ^ A 2 r 2 + (r 2 - m){r 2 - rj) ) ' 

where according to our conventions A is always positive. This expression has a zero at 
infinity, corresponding to a maximum. It has an additional zero at r*, necessarily then a 
minimum, implying = u_ (r* ) . This discussion applies as well to non equal-time geodesies 
in AdS3 with r* > A. 



Substituting in (4.3) and ( 3.16 ) leads to 



1 , 

t Q - — log 
2r h 



m-\{r*+r h f -A 2 )\ , (4.8) 



- 13 - 



Figure 2: Example of a non equal-time BTZ geodesic completely contained in wedge I, i.e. with 



e = e=l and C ( |3.17| ) positive. AdS3 geodesies with r* > A verify the same pattern. 



The value of the constant to is set by the boundary conditions (t4.5|) 

to = ^ + ^log(Cm- 2 ), (4.9) 

and we then have 

h+t 2 , 1 . (r« - r h ) 2 - A 2 

v * = — 7, ^ ~a — 1o St ^ J9- ( 4 -10) 

2 4r h 6 (r* + r h ) 2 - A 2 



Inverting the relation between I and the integrals of motion A and ( 3.21 ), we observe 
that the second term on the rhs equals —1/2. Hence we finally obtain 

l = h + t 2 -2v*, (4.11) 

which generalizes (|4.6|) to geodesies with arbitrary spacelike separated endpoints. For an 
infalling shell with trajectory v = 0, representing a sudden perturbation of the dual CFT 
at i = 0, we reproduce the same threshold for thermal behavior as that found in the study 
of quantum quenches [25,26]. 

4.2 Non-equilibrium geodesies 

We turn now to the construction of geodesies that intersect the infalling shell. The matching 
conditions for geodesies across the shell have been derived in [15]. Let us briefly reobtain 
them here for completeness. Using the affine parameterization, the geodesic equations are 

x = £, (4.12) 

v = -rv 2 + ^| , (4.13) 

r = 2rrv + (r 2 — mO(v))v — —m5(v)v 2 , (4-14) 

where Q(v) is the step function. In the AdS and BTZ backgrounds the geodesic equations 
contained two integration constants, r* and A, associated to invariance under x and t 
translations. The collapse geometry only preserves invariance under translations in x, 
hence above appears only the associated conserved quantity rv The breakdown of time 
translations implies that the energy parameters in the AdS and BTZ segments of geodesies 
that intersect the shell will differ. 
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Let us call A c the value of the affine parameter at an intersection point between the 
geodesic and the infalling shell. We will denote with a subindex "out" the segment that 
lives outside the shell, perceiving a BTZ geometry, and use a subindex "in" for the segment 
inside the shell, feeling AdS3. Continuity requires that 

Xm(K) = x out (X c ) , fin(A c ) = v out (X c ) , n a {X c ) = r out (A c ) . (4.15) 
Then (|4.12| ) immediately implies 

iin(A c ) = Xout(A c ) . (4.16) 

The behavior of v and r across the shell can be obtained by integrating the rhs of 
the geodesic equations between A + > A c > A_ . Unless the integrand contains at least delta 
function singularity, letting A± — > A c will lead to a vanishing result which would imply the 
continuity of the associated first derivative. Notice that a delta function can not appear in 



the first derivative of a continuous function. Applying this reasoning to (4.13) we obtain 



Vin(Xc) = Wout(A c ). (4.17) 

However the first derivative of the radial coordinate acquires a discontinuity due to the last 
term on ( [4.14 ) 



r ou t(A c ) - rin(A c ) = ~mv(X c ) . (4.18) 




(c) (d) 

Figure 3: Examples of different geodesies according to our classification in terms of d; n and d ut- 
a) Geodesic with t\ > 0, t 2 < and d ln = d out = 1; Geodesies with if > and b) d- m — d out = 1; c) 
d in = l and d out = -l; d) d ia = d out = -l. 

Using the matching conditions above, the construction of geodesies that intersect the 
shell can be done in a straightforward although laborious way. We find it convenient to 
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parameterize the geodesies in terms of r*, r c = r(A c ), A m , which denotes the value of the 
energy constant for the geodesic segment inside the shell, and a sign d m describing in which 
branch the inner geodesic segment reaches the shell at r c . The geodesies will intersect the 
shell at one or two points depending on whether t% is positive or negative, i.e. after or 
before the collapse of the shell starts. In case there are two intersection points, they will 
be in general located at different values of the radial coordinate. We choose r c to be the 
smaller one and we will call f c the second, larger value in case it exits. As noticed in the 
previous subsection, when Ai n ^ the minimal value of the radial coordinate reached by 
the inner geodesic segment, r m ( 3.12j ), and the point at which the minimal value of the null 



coordinate happens, r*, do not coincide. Due to this, when r c < r* the intersection point 
can happen in the branch — of the inner segment, corresponding to d\ n = — 1, as shown 
in Fig ||d. The previous data completely determines the shape of the geodesic, but not 
its position in the x direction. Since the system is invariant under translations in x, this 
information does not affect any relevant physical quantity and we will ignore it. 

The only restriction we will assume in the following is r* > A- in . For AdS3 this condition 
corresponds to geodesies with spacelike separated endpoints. It turns out that this choice 
is able to provide all geodesies with spacelike separated endpoints also in the collapsing 
shell geometry. For geodesies whose endpoints tend to be lightlike separated, this fact is 
easily proved. Indeed we have shown that in the lightlike limit geodesies remain tangent 
to the boundary, and close to the boundary the BTZ and AdS geometries coincide. 



From the matching condition ( 4.18|) we obtain the value of the energy parameter in 
the BTZ geodesic segment intersecting the shell at r c 

\A in (2rl - m) - d- m mjA 2 n + r 2 -r 2 
~2r~ 2 



^out = s-g-* • (4.19) 



m 



and the sign d out encoding its associated branch 
d out = e sign 



A 2 n + r c ~ r l + dinA in {m - 2r c 2 ) 



(m - 2r 2 )J A 2 Q + r 2 - r 2 + d in A i: 



rn 



(4.20) 



where e is given by ( p,18| ). Contrary to the thermal case, when the entire geometry is 
described by a BTZ black hole, we need now to allow both positive and negative values for 



the constant C (|3.17|) in order to construct the outer geodesic segment. 

Not all possible data for the inner geodesic lead to geodesic in the collapsing shell 
geometry that are relevant for the calculation of two-point functions. Indeed, some of 
them produce geodesies whose outer BTZ segment falls into the black hole singularity, 
instead of returning to the boundary. In Fig.^ and Fig.|] we have plotted all types of 
shapes that the (real part) of the functions v±(r) in the BTZ background can have 4 . They 
are easily understood using the generalization of ( f4.7| ) to an arbitrary geodesic 

dr r z -m \ y' A 2 r 2 + (r 2 - m)(r 2 - r 2 ) I 



4 The position in the v coordinate is not relevant in this case since it can be shifted without altering the 
shape of the functions due to the invariance under time translations. 
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Figure 4: Examples of BTZ geodesies with a) e = l and C<0; e=— 1 and b) C<0; c) C>0 and 
e = — 1; d) C > and e = 1 , in this case r* < r m . 



When e = l (—1) the branch — (+) has a minimum at r*, having a positive derivative for 
r>r* and negative otherwise. When e = — 1 the branch — always diverges to minus infinity 
at the horizon, having a positive derivative above it and negative below. When e = 1 the 
same pattern for the derivatives holds in the branch +, but the horizon is only reached for 
C < 0. An example of geodesic which should be discarded is in Fig.Qa when d ou t = 1 and 
r c <r^, or when d out = —1 and r c < r*. 

When d out = 1 the only way to avoid that the outer geodesic segment ends at the 
singularity is that v' is positive at the merging point. If d out = — 1 and v' at r c is negative, 
the outer geodesic segment must change into the branch + before reaching the boundary. 
This can only happen in the situations exemplified in Fig.^ and Fig.[|d. Let us define an 
index rj which equals minus one if there is change of branches in the outer segment and 
one otherwise. We then have 




For example, geodesies in Fig.]^ and |3|d have rj = — 1. In Section 6 we will see what is 
the restriction on the data of the inner geodesic segment that ensures geodesies in the 
collapsing shell geometry not to fall into the singularity. 

We will call df n and d^ ut the signs describing the matching branches for the inner and 



outer x(r) functions. The choice (4.4) for the geodesic endpoints sets df n = di n . Combining 



the matching equations ( f!.16| ) and Q4.17 ) we have 

Ct = ^- (4-23) 
Finally, a second intersection point of the geodesic with the shell will exist at 

2 

Tc r * (4.24) 



"2 - 2A in (A in + d in JAf n + r c 2 



r. 



2\ 



whenever the denominator of this expression is positive 5 . For those cases where there are 
two intersection points with the infalling shell, the two separate outer geodesic segments 



5 For r% > 2A? n this is always the case when di n = — 1, while it requires ro < r^/2A- m for di n = 1. When 
instead r\ < 2A? n , the denominator is positive for ro > ri/2Ai n if d\ n — — 1, while it is always negative for 
d in = l. 
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will have in general different energy parameters A out and A out . Also d ou t and d ou t might 
differ. The matching conditions ( [4.16| ) and ( |4.17 ) for the second intersection point lead to 



out 



(4.25) 



We have now all the data necessary to construct the complete geodesic in terms of the 
data of its inner segment. The explicit expression can be found in Appendix A. The proper 
length of the geodesic is f dX = f dr/r, substituting (3.6) with the constants appropriate 
for each segment. The length of the inner segment is given by 



L in = din log [rj (r c + \Jr 2 c - r^)] + log \r n } (r M + J r 



M 



»)]■ 



(4.26) 



When the geodesies intersects twice the inf ailing shell, rj^ = f c . Otherwise tm should be 
replaced by a large but finite value r&, as explained in previous sections. Inverse powers of 
rf, are then to be neglected in the second term of the rhs. The length of the outer geodesic 
segment that starts at r c is 

L out = log 2 r b + ~ log C- 1 (r 2 - A 2 out +m- 2r 2 + 2^ + mr* - r^r* - A 2 out + m)) . 

(4.27) 

In the case that the geodesic intersects twice with the shell, there is an analogous expression 
for the length L out of the outer geodesic segment starting at f c . 




^reg 
10 - 




1 2 3 4 5 6 

(b) 

Figure 5: For m = 1 regularized geodesic length a) as a function of / for fixed At = 3 and 
t\ =3.1, 4.1, 5.1 from bottom up in solid lines {t% > 0), t\ = 0.9, 1.9, 2.9 in dashed lines (£2 < 0) and 
empty AdS geodesic in dotted line; b) as a function of t\ for fixed At = 1 and I = 2, 4, 6, 8, 10 from 
bottom up respectively. 



In Fig.|]a we have plotted the regularized geodesic proper length ( |3.2j ) as a function 
of I for fixed A£ = ti— £2- It diverges when the geodesies approach the lightlike limit. In 
dashed lines we have distinguished geodesies with negative £2 , those associated to two-point 
functions with one time after and one before the perturbation acts. When £1 < the result 
is that of empty AdS3. In solid lines we have plotted the length of geodesic with both 
endpoints after the perturbation. 

As we have shown in subsection 4.1, geodesies with both £, positive and Z<£i+£2 lie 
completely outside the shell. Perceiving a BTZ geometry they reproduce the results at 
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thermal equilibrium. Namely for I bigger than At and such that the associated geodesies 
can reach the black hole horizon 6 , the geodesic length grows linearly 

L reg (Z;ii,t 2 ) = 21og|- + ^-. (4.28) 

A linear regime for L reg can be clearly seen in Fig ||a. Geodesies with I > t±+t2 intersect the 
shell and lead to a result that deviates from thermal. For asymptotically large intervals the 
geodesic length returns to the same logarithmic behavior as in AdS3, plus a time-dependent 
shift 

L TCg (l;t 1 ,t 2 ) = 21ogZ + s(ti,t 2 ). (4.29) 

Plots of the geodesic length as a function of t\ for fixed I and At are shown in Fig.|5]b. 
It saturates to its thermal value at t\ = (I + At)/2, when the associated geodesic stops 
reaching the infalling shell. 

Although our results are is a complicated function of the geodesic data 

which can not be explicitly expressed in terms of the physical variables t\, ti and I. Explicit 
results can be however derived for large intervals or late times. The next sections are 
devoted to analyze the evolution of two-point functions in these two limits, from which we 
will derive important physical information. 



5. Two-point functions for asymptotically large separations 

When I is very large, the associated geodesic will reach deeply inside the shell. In particular, 
for I— >oo the geodesic will get closer and closer to r = 0. We want to take the large I limit 
while keeping t\ and t% finite. Since the major part of the geodesic will be contained inside 



the shell, in the AdS part of the geometry, relation (3.11) implies 



2 , 2a 
— A- ~ — 

v m - 1 



A in ~ -fr , (5.1) 



for some finite constant a, and necessarily d m = 1. The parameters of the outer piece of 
the geodesic are obtained from ( 4.19Q and ( |1.20| ) 



7T& / 9 

Aout = 7T~ » "out = sign(2r c - m) . (5.2) 
2r c 

If the first crossing happens at the value of the radial coordinate r c , the second should take 
place at f c = r c /(l — ar c ). Hence there will be a second intersection point only if a<l/r c , 
in which case expressions analogous to ( |5.2| ) hold. 

The length of the AdS geodesic segment, given by ( [4.26 ) , is 



L in = 2 log I + log r c + log r M + O^ 1 ) , (5.3) 

where as before tm equals f c if the geodesic intersects twice with the shell and r& otherwise. 
Here / strictly represents the space separation between the endpoints of the inner geodesic, 
prolonged to the boundary in the absence of the shell. However the difference between 



"This distance is of order f3/2 for equal-time geodesies. 



- 19 - 



this separation and that calculated in the presence of the shell, having into account the 
matching conditions, is a finite number. Thus (|5.3| ) applies equally with I being the space 
separation between the geodesic endpoints in the thin shell geometry. The length ( |4.27| ) of 
the outer geodesic segment starting at r c and its endpoint time t±, given by ( A.5| ), are 



L out = log2r fc -log^^, t 1 = llog^±^, (5.4) 

with analogous expressions holding for the second crossing at f c . We have used the fact 
that r c >rh/2 for r„,~0, which will be explained in the next section. The relation between 
ti and the associated intersection point can be easily inverted. Adding the length of the 
several geodesic segments and subtracting 21ogrf,, we finally obtain 

L reg = 2 log I + s(tx) + s(t 2 ) + Oil' 1 ) , (5.5) 

where 

m= f21og[cosh(7rf//3)], t>0 
W \0, t < ' { ' ' 

Comparing with ( |4,29| ), we have s(t\, t 2 ) = s(t\) + s(t 2 )- This precisely matches the numer- 
ical results of [9]. 

Substituting ([)!]), we obtain the asymptotic value of the two-point function 



<0(Z,ti)0(O, h)) ~ — ^ . (5.7) 

(I cosh(vrti//3) cosh(vrt 2 //3)) 

Its dependence on / reproduces that of the two-point function on the CFT vacuum for 
large spatial distances. On the contrary, an exponential suppression with / was found in 
the same regime after a quantum quench, see (|2.8| ) [25,26]. Both quantum quenches and 
our holographic model generate at i = a sea of excitations with non-trivial entanglement 
properties. From then on coherent excitations start separating from each other. However 
if we consider very large for any finite ti they will have not separated enough to erase the 
effect of the initial entanglement pattern on the two-point functions. Hence the different 
behavior of the two-point functions in both dynamical situations should be traced back 
to the different initial states that trigger the evolution. In the holographic set up, the 
perturbation modeled by the shell does not destroy the long range correlations proper 
of the CFT vacuum [9]. However the ground state of a massive hamiltonian was taken 
as initial state for the quenches considered in [25,26]. In an state with a mass gap, the 
exponential damping of correlators with the distance is expected. 

Another important feature of ( |5.7| ) is its increasing suppression as a function of time 
due to the contribution of the s(ti). This function has a linear behavior for times t>/3 

.s(t)^-21og2 + -^, (5.8) 

as noted in [9] . Let us denote i the threshold time from which on the previous approxi- 
mation is valid. The time dependence of the two-point function (|5.7l ) for ti > t reduces to 
e -27rA(ti-w 2 )/^_ Identifying (3 = 4tq, this is in agreement with the quantum quench result. 
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6. Steady evolution 



For finite separations l>t\ + t2, the I dependence of two-point functions after a quantum 
quench turns out to enter in terms of the effective distance l—t\—ti, as seen in Q2.8| ). We 
will show in this section that the same pattern applies to the holographic model for times 
after the threshold i. Namely, the power law behavior of two-point functions with / that 
we found in the previous section for very large distances, gets promoted for finite l>t\+t2 
and U >i to a function of l—t\—t<2. The evolution of the entanglement entropy in the same 
holographic set up verifies an analogous property. It was shown numerically in [9] that the 
entanglement entropy of a single interval of size i, for l>2t and t after the same threshold 
time i, is a function of the combination I — It. 

In order to arrive to the above conclusion, we first need to recall some further properties 
of the geodesies. As mentioned in Section 4, not all the initial data r*, r c , A- m and d\ n lead to 
geodesies relevant for the computation of two-point functions. Indeed some of them produce 
geodesies that fall into the black hole singularity instead of returning to the boundary. The 
set of initial data that avoid this problem is 



r* > r h + A h 



r c £ [r m ,oo] , d in = 1 



r c £ [r m ,oo] , d in = 1 ( 6 - 1 ) 
r c £ [r m ,R] , d in = -1 



r* G [f*,r h + A in ] 
r* 6 (A in , f *] : r c £ [R, oo] , d in = 1 . 



where 



r h + 2A in - 2r* + y/(r h + 2A in ) 2 + 4r,(r h - rv)) 

^ = ^ 4(r, + A n -r,) • ^ 

and f* is the only value of for fixed Ai n at which R = r m , with r m the minimal radius of 
spacelike AdS3 geodesies ( p,12| ). The functional dependence of r* on A- m is given in (B.J). 



For any other r#, i? > r m . Consistently at r c = rh + A[ n , R = r*. See Appendix B for a 
discussion of the properties of the functions R and f*. 

We have derived the bound on the intersection point by inspection over the param- 
eter space of the geodesies. The casuistic for the construction of generic geodesies is so 
broad that unfortunately we do not have a closed proof for fl6.2|) . An illustrative exam- 
ple for the necessity of introducing R is provided by r* = A m = 0. From Q5.2D we know 
that when r c <r^/V2 the inner geodesic merges at the shell with the branch — of a BTZ 
geodesic. Using ( f4.22|) we have that for these values 77 = — 1. For r c = r/ l /2, the minimal 
radial value ( 3.19| ) at which the branches ± of the BTZ geodesic join happens at the sin- 



gularity. For smaller values of r c , e becomes negative implying the outer geodesic will have 
two disconnected pieces each of them reaching to the singularity. Hence r c must be bigger 
than R = rh/2 in order to have an smooth outer geodesic that returns to the boundary. 

A very important property of ( |6.2| ) is that as r c ,f c ^-R, the functions tj, I and L reg 
turn out to diverge. Let us define 5r c = \r c — R\/R and 5f c = \f c — R\/R. In order that both 
r c and f c tend to R, definition (4.24) implies that A in must be small. In particular, A- in 



-21- 



must be of the order of 5r c and 6r c . We show in Appendix B that in this limit 

2r h (r h - n)(B% - r h r*) 



2^ l0g ^+2^ l0g 



i log5r " c+ i: log 



r*fo(r h - Ro) 2 
2r h (r h - n)(R% - r h r*) 



.— lo g (8r c 5f c ) + — log 



r*fo(fh - Ro) 2 
2r h r*(Rl - r h r*) 



+ 



4(r/, - r*)R 



(6.3) 



L reg = -- \og(5r c 5r c ) + log 



2(^§ - r h i 



+ 2 log 



Ro + VR^ 



2 rl 



r* 



forhr*(r h - r*) 

where i?o is the value of i? (|6.2D for yl; n = 0, and /o a function of r* whose expression will 
not be relevant for us. These relations are valid up to terms of order 5r c , 5f c and A{ n . 
Combining the previous equalities, we obtain 



, , , 2 (r h - R )n 4(r h - r*)R 

I - t x - t 2 = — log - — H 

rh {rh -r*)Ro mr* 



(6.4) 



The importance lesson here is that the rhs is only function of r*. From (6.1) we know that 
R only plays a role for r* <rh+Ai n . Since we are taking the inner energy parameter to be 
zero, we must consider r* in the interval [0, rh]- Consistently, this range covers univocally 
values of the lhs from just the threshold for thermal behaviour at r*=r/j where we have 

l-h-t 2 = 0, (6.5) 

to infinite spatial separation for the geodesic endpoints as r*— )-0, when (|6.4[) leads to 



I 



as previously derived in (|5.1|). 

Using (|6.3D , we obtain for the geodesic length 



(6.6) 



Lreg = r h (h + t 2 )-2 log r h + 2 log — + 2 log ■ 



(6.7) 



rh ~ r* r* 

At r* = rh, the threshold for thermal behavior, this expression reproduces the extensive 
dependence (4.28]) characteristic of thermal equilibrium. At r*— >0 we have 



L 



reg 



r ft (ii + i 2 )-41o g 2 + 21ogZ. 



(6.8) 



Agreement with the result ( |5.5| ) for the geodesic length at large space separations requires 
ti > i. Extracting the values of ij for which ( |6.3| ) applies is involved. However using 
r/j = 2ir/ft, a rough analysis indicates that they are consistent with the bound set by i. 
Hence, inverting ( |6.4j) to obtain r* as a function of l—ti—t 2 , we derive the announced result 

2vr 



L 



rcg 



ft 



(h + t 2 ) + F(i-ti-t 2 ), 



(6.9) 
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for ti>i and separations I > t\ + t2- The validity of the above expression can be further 
checked in Fig. ^a, where we have plotted the geodesic length as a function of l—t\—t2 for 
several values of t\ keeping At fixed. We observe that as t\ increases the curves tend to 
coincide up to a shift, as implied by (plj|). 

We can identify two sources that drive the evolution of an out of equilibrium system 
towards equilibration: the interaction between its components and their propagation. In 
the sea of excitations sourced by the initial perturbation, quantum correlations will be 
stronger between excitations generated at close points. Propagation will tend to increase 
the separation between entangled excitations with time. On the other hand interaction will 
lead to the redistribution of energy, and other conserved charges in case they are present, 
among the different momentum modes. This is crucial in order that the state towards which 
the system evolves is that of thermal equilibrium. The separation of entangled excitations 
with time is a process constrained by causality, implying that equilibration can only be 
achieved on regions of growing but finite size. On the contrary, the redistribution of energy 
among modes can be expected to require a finite time for its nearly completion. 

In Section 2 we have reviewed the simple model proposed in [23] to account for the 
evolution pattern of the entanglement entropy and two-point functions after a quantum 
quench. It explained the dependence of two-point functions on the effective distance l—t\—t2 
for separations I > 1 1 H- ^2 • The main hypothesis of this model was that the excitations 
generated by the perturbation move as free particles according to their group velocities, 
which for a 2-dimensional CFT is always given by the speed of light. Therefore it seems 
consistent to assume that when this model is applicable, all processes in the evolving plasma 
related with interaction have reached equilibrium. In this sense, notice that the general 
results for the evolution after a quantum quench derived in [23, 25, 26] are only valid for 
times much larger than the inverse effective temperature. As it was proposed in [9] based 
in the study of the entanglement entropy, this suggests to interpret t in our holographic 
model as the thermalization time for occupation numbers. Here we have confirmed the 
same conclusion from the evolution of general two-point functions. 

It is interesting to cross check our arguments with the simple example of a quench from 
a massive to massless free boson, studied in [26]. In this case the occupation numbers are 
conserved along the evolution. They do not coincide with a thermal distribution, except for 
momenta smaller than the effective temperature T e ff = m/4, where m is the mass before the 
quench. This provides a model where equilibration is not synonymous of thermalization. 
Since the occupation numbers are conserved we should have if ree = 0. An exact treatment 
is possible for all times and distances, and the following relation turns out to be fulfilled 

(0(t, l)0(t, 0)> = e -A(M(t)+F(l-2t)) ( 610 ) 

for any time after the quench and distances / > /3 e // 7 . Notice that also ( |6.9[ ) requires a 
7 For the free boson quench we have 

§(i) = -i/(2mi), F(x) = jmx-± + ±f(mx), (6.11) 

where /(z) = 1 + §G?3 j () 3 ( 2 \ with G the Meijer G-function. 
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(a) (b) 

Figure 6: For m = 1, a) L* cg = L rcg — L rcg (l = £1+^2) as a function of I — 1\— 1% for t\ — 
1.5,2.5,3.5,4.5,5.5 from top down and At = 1. The last two lines practically coincide with the 
function F in ( |6.9D (up to a constant), b) L rcg as a function of I for At = and £ = 0, 1.5, 3, 4.5, 6, 7.5. 
We have signaled approximately the times to and t. 

restriction on the distances where it applies. The function s for the free quench ( |6.11| ) has 
a similar structure to that of the holographic model. It vanishes at t = and after some 
time it enters a regime of linear growth. The main difference with the holographic model 
is that the time at which s behaves linearly does not coincide with if ree . This example 
suggests that the important property which could identify the equilibration of occupation 
numbers in the models we are considering, is when the dependence on the scale I enters 
through the effective distance I— 2t or l—t\ — t2. 

In the free boson model the time at which s starts growing linearly coincides with 
that at which L reg begins to show an extensive behavior for some range of I. Let us call 
the time defined by this second property tQ. In the holographic model io~/3/4. It is the 
time at which equal-time geodesies anchored on the boundary start reaching the black hole 
horizon and produce a result clearly identifiable as thermal. It can be associated through 
propagation at the speed of light with the smallest scales able to feel the thermal bath. 
The time i is clearly bigger than to in the holographic model (see FigJ^b). The fact that 
the latter can be explained from propagation provides an additional argument in favor to 
relate the former with interaction. 

7. Discussion 

In this paper we have studied the evolution of two-point functions in a 2-dimensional CFT 
using holographic techniques. In our model the CFT starts in its vacuum state and at 
t = undergoes an instantaneous perturbation that brings the system out of equilibrium. 
The dual gravitational background that represents this process is that of the collapse of an 
infinitely thin shell of null dust in a spacetime with negative cosmological constant. It has 
been argued in [9,24] that this set up represent a unitary process in the field theory, implying 
that information is not lost at the microscopical level. However very general arguments [28] 
show that the system will evolve towards an state that at the macroscopical level can be 
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described as thermal equilibrium. This is reflected hole-graphically in the formation of a 
black hole out of the gravitational collapse of the null dust shell. 

An important question is how is the energy is distributed among modes of the initial 
excitations created by the perturbation. Both the Vaidya metrics and their thin shell limit 
asymptote at large radius to AdS,3, differing from it only for values of the radial coordinate 
set by the mass function of the black hole being formed. Since this is the only scale in 
the problem, it suggests that excitations are produced around energies of order \J m(t) 
at a rate of the order dy/m(t) /dt. Hence UV excitations would be much suppressed in 
the initial state [9] 8 . This is supported by the fact that both entanglement entropy and 
two-point functions never deviate from their vacuum results in the UV limit. Recently 
a holographic set up which allows flexibility on the initial conditions for the evolution 
has been proposed [16]. It would be very interesting to study the dynamics of non-local 
observables in these cases. 

Although the previous reasoning indicates that the initial excitations are probably not 
produced very far from their equilibrium distribution in our model, a deviation from it is of 
course to be expected. Based on the holographic analysis, a time t after which the energy 
redistribution among modes according to thermal equilibrium should be nearly completed 
can be proposed [9]. One of the main result of this paper is that the same time i derived 
from the analysis of the entanglement entropy, applies to generic two-point functions. After 
t and for l>t±+t2, the dependence on the spatial separation of two-point functions enters 
only through the effective combination l—t\—t2 in accordance with the late time evolution 
after a quantum quench [25,26]. We have argued that this is the relevant property to 
identify the equilibration time for occupation numbers in these 2-dimensional systems. 

Both in quantum quenches and in our holographic model the initial sea of excitations 
has a non-trivial entanglement pattern, where the main difference between both cases is the 
presence of long range entanglement in the latter. In order that observables over an scale I 
can reproduce a thermal result, it is necessary that entangled excitations have separated at 
least a distance I. Causality then implies that thermalization can not be achieved globally 
over the system. For 2-dimensional CFT's, all excitations propagate at the speed of light. 
This gives rise to a sharp cut between the size of regions over which thermalization is 
effective and those in which it is not. A second central result of this paper is to reproduce 
the condition for thermal behavior of general two-point functions, I < <i+t2 with > 0, 
first derived in the study of quantum quenches [25,26]. 

When the dynamics of an out of equilibrium system is unitary, a very interesting 
question is how to define a coarse grained entropy that gives a good description of its 
evolution towards macroscopical equilibration. In [24] it was proposed that an interesting 
definition of a coarse grained entropy for compact systems could be given in terms of the 
entanglement entropy of a region A, where A covers half of the total system, through 



with the entanglement entropy in the state before the evolution starts. A natural 
extension of this relation to the non-compact set up we are considering is to substitute 
8 This is also the case in the free boson quench [26] 




(7.1) 
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the region A by the complete non-compact system. Both for the 2-dimensional quantum 
quenches studied in [23] and the holographic model based on the infalling shell of null 
dust [9], the entanglement entropy of a single asymptotically large interval satisfied 

S oo (t) = ^s(t) + S oo . (7.2) 



In the infinitely thin shell limit of the holographic model s(t) is given by (|5.6| ). For the 
critical quantum quenches this function is only known for late times, where it coincides 
with fl5.8|) up to possibly a numerical factor. Therefore the generalization of ( |7.1|) to the 
non-compact case would be 

Seff(i) = • (7.3) 

This definition implies that for any finite time S e g is finite. Hence no finite entropy den- 
sity is generated, in agreement with the fact that equilibration is only effective over finite 
regions. Notice that this rules out the apparent horizon, in the theories with a holographic 
dual, as the place where the field theory entropy resides. Moreover in the thin shell model 
the apparent horizon, which is located at r = y/m@(v), would predict an instantaneous 
thermalization for the entropy. On the contrary the function s shows an smooth interpola- 
tion between zero and the thermal density multiplied by the size of the largest thermalized 
region for late times. 

As we mentioned in Section 6, the function s for the free boson quench ( |6.11|) has 
a similar structure to that of the holographic model, which describes thermalization in a 
strongly coupled plasma. The main difference between them affects the time at which they 
start to behave linearly. It would be very interesting to better understand the information 



about an equilibration process that can be extracted from (7.1) and 



Figure 7: Geodesies whose outer segments are of the type that would connect the two boundaries 
of the extended BTZ Penrose diagram. They enter well behind the apparent horizon. We have 
chosen t\ = 6 and a) At = 0.4 and b) At — 1.4. In a dashed line we have plotted the apparent 
horizon 

In the Schwinger-Keldysh or real time formalism, the thermal entropy can be under- 
stood as entanglement entropy in the pure state (3.25) defined on two copies of the system, 



when tracing over one of them. In [38] the two copies of the system were put in correspon- 
dence with the two asymptotic regions of the extended Penrose diagram of the AdS eternal 
black hole. In this extended diagram there are geodesies that connect the two boundaries, 
reaching behind the horizon without falling into the singularity. Namely, those for which 
the constant C (|3.17|) is positive and the signs ( p.18 ) satisfy e = — 1 and e = 1 (see Fig.||). 



The function s in the holographic model is derived from the length of geodesies whose 
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endpoints have a large spatial separation, corresponding to small value of the integral of 
motion r*. Relations (5.1) and ( |5.2| ), together with (5.4), imply that from quite early on 



in the evolution the BTZ segments of these long geodesies are precisely of the mentioned 
type (see Figj^for examples). This offers a very consistent picture. The coarse grained en- 
tropy that describes the approach to thermalization requires information from behind the 
apparent horizon. Moreover this information is extracted from geodesies that in the static 
black hole would join the two boundaries, strongly reminding to the real time formalism. 
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A. Geodesic trajectories 

Below can be found the explicit expressions for the different geodesic segments in the 
collapsing shell geometry. When dj n = l, the inner geodesic segment is given by 



r € [r m ,r c ] 



r G k m ,r M J 



v(r) 


= t +t AdS (r) 


- 1/r 


x(r) 


= x AdS (r) 




v(r) 


= t -t AdS (r) 


- 1/r 


x(r) 


_ _ x AdS^ 





(A.l) 



where r m given by ( 3.12] ), and rj\/ = r c if there is a second crossing of the geodesic with 



the shell and infinity otherwise. The functions t and x are given in ( p.lOj ). When 
instead di n = — 1, we have 

v(r) = t — t AdS (r) - 1/r 
x(r) = —x AdS (r) 

The integration constant to has to be chosen such that i>d in (r c ) = 0, and for simplicity we 
have set xq = since it does not affect any relevant quantity. 

The outer geodesic segment starting at the intersection point r c has the following 
expression. When ?j=l we have 

v(r) =h- t dout (oo) + v dout (r) 
x(r) = x + (r) — x + (r c ) + di n x Adb (r c ) 

with td and Xd given by ( ^.16 ). For 77 = — 1 we must have d out = — 1 and then 



re [r c , r 




(1)1 J v ( r ) = h - t+{oo) + v-(r) 



x_(r) — x_(r c ) + d- m x AdS (r c ) 

(A.4) 



(l) 1 I v(r) = ti - t+(oo) + v + (r) 



r € [r^/,00] : 



x(r) = x + {r) — x_(r c ) + d- m x AdS (r c ) 
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with the minimal value of the radial coordinate defined in ( 3.1S| ). The value of the 



time coordinate at the endpoint of this geodesic segment, ti, is given by 

h = t v d out (oo) - v dout (r c ) . (A.5) 
The outer geodesic segment starting at the intersection point f c , when 77 = 1 is given 



by 



r G [f c , 00] : J d{r) = h ~ (00) + (r) , (A.6) 
I x(r) = X-(r) — x_(r c ) — x AdS (r c ) ' 



where the bar over the functions indicates that the energy parameter should be taken to 
be A out instead of A out . For fj = — 1 we have 



(2) 1 J v(r) = t 2 - t + (oo) + v-(r) 

(A.7) 



r G [r c ,r m j . % _ _ _ AdSf- , 

x(rj = x_|_(rj — x + (r c j — x (r c j 



r G [rV , 00] : 

with 



( 2 ) ~0 . J u ( r ) = *2 - *+(oo) + v + (r) 

x(r) = X-(r) — x+{f c ) — x AdS (f c ) 



t^kdoJ^-v^Jfc). (A.8) 

From the previous expressions we can also read immediately I, the spatial distance 
between the geodesies endpoints. If there are two intersection points with the shell we 
have 

I = x+(oo) - x v (r c ) + d in x AdS (r c ) + x AdS {f c ) + x+(oo) - Xff (f c ) , (A.9) 
where we have used that x_^(f c ) — x_(oo) = x + (oo) — Xfj(f c ). Otherwise 

I = x+(oo) - x v (r c ) + d in x AdS (r c ) + x AdS (oo) . (A.10) 
B. Geodesies in the steady regime 



We will start deriving several properties of the function R ( |6.2| ), which bounds the allowed 
values of the intersection point r c with the shell for geodesies that do not fall into the black 
hole singularity. The function R depends on r* and A m . Its derivative with respect to r* 
is always positive, and thus it is an increasing function of this variable for fixed A- in . By 



substituting in (6.2), it is easy to verify that 

r* = r h ,r h + A in => R = . (B.l) 
These are the only points where this property is verified. Also easy to check is that 



dR 

dr* 



r h r h + A in 



which implies 

r h < r* < r h + A in =>- R < r* , < r h =>• R > . (B.3) 
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Between and r^ + A^ there is a single f* at which R = r m , with r m the minimal radial 
coordinate of spacelike AdS3 geodesies. The value of r c is given by 



r h +A in + 



2 2 / 3 (r 2 h + 2A in r h +AAf n ) 



(AtI - 15 A in r 2 - 24A 2 Q r h - 32 Af n +3^-84^-13^-164? 



^3 
in /i 



1/3 



+ 



(4j - 154n?j - 24Af n r fe - 32^+3^-84^ - 13^ - WA* 

22/3 



1/3 



, (B.4) 



The fact that r* > r m together with ( |B.1| ) , imply that R is bigger than r m for any other 
r*. All these properties are resumed in Fig.||a. In Fig.||b we have plotted f*— 4n- This 
is a decreasing function of the energy parameter that coincides with the horizon radius at 
vanishing 4n and tends to for large values of 4n- 





Figure 8: For r^ = l, a) plot of R, r* and r m for 4n = -5; b) plot of r*(4n) — A- ln . 

We will study now the properties of geodesies whose intersection points with the in 
falling shell are close to R. By direct inspection, one can check that the parameter C (3.17 
vanishes at r c = R, implying that 



A 2 out +r 2 h \ R = 2e(R)r,r h . 



(B.5) 



Since A 2 ut ( |4.19| ) is a continuous function of r* over the whole interval [4n, + A- m ], 
i(R) must take a unique value. We can therefore check it at a particular r*, which we 
conveniently choose to be the upper limit of the interval. At r* = rfc + 4n the inner and 
outer energy parameters coincide, 4ut = 4n, and hence e{R) = 1. Substituting back in 



(B.5), when r c —^R we then have 

4>ut = Vh - r*\ + sign(r h - r*)f8r c 



(B-6) 



where 5r c = \ r c — R\/R and / is a function whose explicit expression we will not need. In 
the same limit we can check that 



sign^ - r*J 



(B.7) 
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Repeating the continuity argument above shows that the function whose modulus 
defines ^4 Q ut verifies at r c = R 



A m (2R 2 - m) - d in myjAl + R 2 - r 2 
2iP 



±(r* 



rh) 



(B.8) 



where one sign on the rhs must apply to the complete interval where R is defined. Selecting 
again r Jf = rf l -\-A- m , we easily obtain that the correct sign is plus. Substituting this relation 
in (4.20) we derive after some algebra 

cUt = - sign ((to - 2R 2 )(r h + A in - r*) + mA- m ) . (B.9) 

It can be checked that the zeros as well as the sign of the previous function coincide with 
those of the simpler function R 2 — r^r*, and hence 

d out = -sign(R 2 -r h n). (B.10) 

Finally, using ( 4.22 ) and the properties of R summarized in FigJ^a we deduce 

r ] d ont = l, (B.ll) 

Analogous expressions hold for the second crossing whenever f c ~i?. 

The previous results provide all the necessary ingredients to derive the values of ti, I 
and L reg in the limit that the intersection points with the shell are close to R. We start 
by analyzing the value of the time coordinate at the geodesic endpoint t\, given by (|A.5| ). 
Using several of the results above, we get 



hd out (oo) 



1 



log (|r^ - (r^sign^-r^-Aout) 2 ^ x ) . 



(B.12) 



From ( |B.6| ) we see that the argument of the logarithm vanishes as r c tends to R. However 
Vd out (r c ), also contributing to t\, remains finite in that limit 



1 A\{r h - n )(R 2 



rh{rh ~ Rf 



(B.13) 



As explained in Section ||, in the limit that both intersections points are close to R, the 
energy parameter of the inner geodesic segment A- m must be small and of the order of 5r c . 
We are interested in obtaining the value of t\ up to terms that vanish in the limit r c , f c —}R. 
Therefore it makes sense to set Ai n = in the finite contributions to t\. Combining the 
previous two pieces with this further approximation, we obtain 



1 1 r 1 , 

-— log 5r c + — log 
2r h 2r h 



1r h (r h - r*)(i?Q - r h r*) 



(B.14) 



where Rq and fo are the values of R and / for ^4; n = 0. The same expression applies for 
replacing Sr c by Sf c . 
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Several pieces contribute to the spatial separation I of the geodesic endpoints, given 
by ( |A,9| ). Prom the outer BTZ geodesic segment starting at r c we have: 



1 



x+ (oo) = -— log 
2r h 



\Al 



- {rh-r*) 2 \m x ) , 

which is again divergent as r c tends to R, and x^frc), which remains finite 



(B.15) 



(B.16) 



There are analogous pieces from the crossing at r c . The geodesic segment inside the shell 
contributes to I as di n x AdS (r c ) + x AdS (f c ), with x AdS given in ( |3.9[ ). This quantity is finite 
when the intersection points tend to R and, as explained above, we can set A\ n = to the 
order we are working. When A{ n = the minimal value of the radial coordinate attained 
by the inner geodesic equals r* and, as seen in Fig.|8|b, r* = in that case. The range of 
allowed data for the geodesies ([O]) implies then that r c can only happen in the branch 
d m = l. Therefore the contribution of the inner geodesic to / is 2x AdS (R$). Collecting all 
these pieces we obtain 



I = -— log(5r c 5r c ) + — log 
2r h 2r h 



2r h n(Rl 



rhr*) 



+ 



4(r fe - r*)Ro 



(B.17) 



The geodesic length gets a divergent contribution from the piece logC in ( 4.27 ) in the 
limit in which the intersection points tend to R, given that in this limit C vanishes. Setting 
further ^4; n = and substituting in ( |4.26 ) and (4.27), we have 



L 



reg 



^\og(5r c 5r c ) + log 



2{Rl 



ThU) 



forhr*(r h - r*) 



+ 2 log 



(B.18) 
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